APS/123-QED 



Formation of the internal structure of solids under severe action 
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On the example of a particular problem, the theory of vacancies, a new form of kinetic equations 
symmetrically incorporation the internal and free energies has been derived. The dynamical nature 
of irreversible phenomena at formation and motion of defects (dislocations) has been analyzed by a 
computer experiment. The obtained particular results are extended into a thermodynamic identity 
involving the law of conservation of energy at interaction with an environment (the 1st law of 
thermodynamics) and the law of energy transformation into internal degree of freedom (relaxation) . 
The identity is compared with the analogous Jarzynski identity. The approach is illustrated by 
simulation of processes during severe plastic deformation, the Rybin kinetic equation for this case 
has been derived. 

PACS numbers: 05.70.Ln; 05.45.Pq 



I. INTRODUCTION 

Thermodynamics of internal state variables has a long 
history. An important idea to take into account internal 
microstructure changes for polyatomic gas by introduc- 
ing an additional variable was proposed by Herzfeld and 
Rice in 1928 yy]. Later, the concept was used by Man- 
delshtam and Leontovich for solving the same problem 
for liquids Supplemented with elements of rational 
mechanics the concept took a completed and closed form 
in the works by Coleman and Gurtin Q. The concept 
was given an alternative development in the works by 
Landau. On the base of newly introduced internal state 
variable in a form of the order parameter a theory of 
phase transitions was constructed by him at the quality 
level 0,|5]. 

Studies of the problem have proposed a large variety of 
kinetic equations describing changes in the internal struc- 
ture under nonequilibrium conditions within the concep- 
tual framework of constitutional theories 0, 0] , contin- 
uum theories of defects jsj, [13, |Tl[ , extended contin- 
uum theories 13| , self-organization theories of dislo- 
cations llq, gradient [TtI . [T^ and non-local the- 
ories J^, mesoscopic theories 2l|, [l^], and so on. 
There is a special need to mention the works of Landau's 
School [1, Q ) which was developing one of the most phys- 
ical lines of the theory of non-equilibrium phenomena. 
In the modern science the line is presented by phase- 
field theories [H [H HE [H, [13, H, HI . Unlike Landau 
theory, the order parameter of phase-field theory is not 
well-defined, but it does not hinder the obtaining results 
which are in a good accordance with the behavior of real 
systems. 

The order parameter introduced by Landau is a type 
of the internal variable and therefore, it should describe 
changes in real structure state of a system. However, it 
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describes the state in terms of structure itself not directly 
but through an immediate indicator connected with sym- 
metry change. It may turn out that in some applications 
it is more convenient to directly connect a system state 
with the change of a specific structural parameter, for 
example, the density of defects. 

Internal phenomena proceeding in a solid represent a 
complex pattern of mutual transformations of the energy 
from each type to other ones. In most cases the physical 
nature of internal phenomena is not defined concretely 
[30| . It is assumed that all of them are of fluctuation ori- 
gin [3ll.[3^. At the same time, the generation and motion 
of structural defects in a solid under a severe external 
action can be also regarded as some internal processes 
characterized by a set of internal variables [33]. These 
channels of energy dissipation still differ from traditional 
viscous channels of dissipation as the relaxation involv- 
ing the defects has the dynamical nature. It doesn't run 
directly by means of transition of the organized energy 
to chaotic thermal one but passes a series of intermediate 
levels and stages. During generation of defects a part of 
organized (elastic) energy is spent directly for the forma- 
tion of structural defects which, along with the thermal 
channel, make parallel channels of dissipation. 

The other part of the energy transforms to thermal 
form not immediately, it is initially radiated in the form 
of low-frequency vibrations and waves (acoustic emission) 
which can be treated as a part of the general dissipation 
phenomenon [H, [H, [s^. Thereupon this part of the 
energy is scattered on the equilibrium thermal vibrations 
(waves) replenishing their energy. So, the energy stored 
in the non-equilibrium subsystem cannot be too much. 

The present paper is devoted to the development of 
the Landau approach. In Sec. II on the example of 
the vacancy theory we propose a new form of kinetic 
equations for defects which is symmetrical with respect 
to free and internal energy utilization. In this section 
invoking molecular dynamics simulation, we select equi- 
librium and nonequilibrium heat subsystems of a solid 
under severe load and introduce a notion of nonequilib- 
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rium temperature. In Sec. Ill the results obtained in the 
previous sections are generahzed in the form of a thermo- 
dynamic identity combining the first and second laws of 
thermodynamics for the described processes. In Sec. IV 
the general results are discussed. Section V is devoted to 
the application of kinetic equations to the description of 
grain refinement during severe plastic deformation. Our 
conclusions are presented in Section VI. 



II. KINETIC EQUATIONS 

To begin the analysis of the problem it is convenient to 
maximally simplify it. Let us split the problem into two 
parts, the first, structural one connected with generation 
of defects and the second, thermal part, and consider 
them separately. Moreover, we consider the structural 
part of the problem on the simplest example of a solid 
with vacancies. 



A. Structural part of the problem 

A classical theory of vacancies was developed by 
Frenkel and Kittel [s^, • In the variant proposed by 
Frenkel, the internal energy of a solid with vacancies is 
presented in the form U = Uq + uyn, where Uq is the 
internal energy of the vacancy-free crystal, uy is the av- 
erage vacancy energy, n - is the number of vacancies. 
Next, the transformation to the configurational free en- 
ergy Fc = U — TSc, where T is the temperature, Sc is 
the configurational entropy, is performed. The configu- 
rational entropy is uniquely expressed through the num- 
ber of vacancies, and, at low vacancy density, it is true 
Sc = kN [InN / n + 1), where k is the Boltzmanns con- 
stant, N is the total number of lattice sites occupied by 
atoms. The equilibrium value of the number of vacancies 
Ue is found from the free energy minimum provided the 
vacancy energy is independent of the number of vacan- 
cies, that is, uv = wyoj where uvq is a constant value for 
a given matter. The free energy minimization procedure 
leads to the equation of state rig = Nexp{—{uvo/kT)). 

With the vacancy energy dependence on vacancy num- 
ber taken into account we consider a quadratic correction 
to the internal energy in the form: 



U = Uo + 



1 2 



(1) 



where uyi is a coefficient. Since the energy needed for 
the formation of a new vacancy, in conditions when in 
a solid there are already a number of other vacancies, is 
less than in the case of the vacancy-free crystal, then the 
negative sign must be used for the correction term. Note 
that expression (1) is true both for equilibrium and non- 
equilibrium states. In this approximation the internal 
energy is a convex function of the defect number hav- 
ing the maximum at point n = Umax, as it is shown in 
Fig. [T^. For the transition to the free energy with the 



Legendre's transformation, it is necessary to subtract the 
bound energy TSc from the internal one. According to 
Kittel, a theory of vacancies can be equivalently formu- 
lated either in terms of the configurational entropy or 
the chemical potential [38]. In this context, the bound 
energy can be written down in two equivalent forms: 



TSc 



Uyn. 



(2) 



Strictly speaking, relationship ^ is just only in the 
equilibrium form, as the equilibrium temperature T is 
contained there. Following the Landau style, we in- 
troduce the generalized thermodynamic force for non- 
equilibrium state in the form: 



Uv 



dU 
dn 



uvo — uvin. 



(3) 



Relationship ([3]) can be treated as an equation of state 
true for both the equilibrium state and non-equilibrium 
cases. Applying, though formally, the Legendre transfor- 
mation with respect to a pair of the conjugate variables 
Uv and n and with the equation of state ([3]) one obtains 
for the free energy: 



Fc = Un 



1 



2uv: 



■{uvo ~ uv^ 



(4) 



Note that an eigen argument of the internal energy is 
the vacancy number (density) , and the eigen argument of 
the configurational free energy is the vacancy energy. In 
this approximation the free energy is a concave function 
with the minimum at point uy — uvoj as it is shown in 
Fig. [Ih. 

We define the equilibrium state from a macrocanonical 
distribution taking it in the form: 



fin) 



C(^exp(- 



kT 



(5) 



The pre-exponential factor describes the combina- 
tional, that is, entropy, part of the distribution func- 
tion connected with degeneration of microstates. The 
exponential factor describes the restrictive part of the 
distribution function connected with the overcoming of 
potential barrier between microstates. The most proba- 
ble state determined from condition df{n)/dn = gives 
the equation of state in the form: 



Uv . 



iVoxpl-^), 



(6) 



which is true only for the equilibrium case. As for the 
equilibrium state the equations of state ([6]) and ([3]) must 
be identical, then substituting Eq. ^ into ([3]) one gets 
a condition for the determination of this state: 



N ,kT 
kTln — = Uvo + (7777 - uvijUe 

He 2.N 



(7) 



With relationship ([7]), it is easy to show that neither 
the maximum of the internal energy nor the minimum of 
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the free energy satisfy this state. The equihbrium state 
is at point n = Ue, where relationships 



UVe 



dU 

dUe 



OF, 
du 



Ve 



(8) 



are vahd. Here the additional subscript e denotes the 
equilibrium value of a variable. 

If a system has deviated from the equilibrium state 
then it should tend to that state with a speed which is 
the higher, as the deviation is larger: 



dn 
'dt 



dU_ 

dn 



UVe), 



duv 
~dt 



dF 



The both variants of the kinetic equations are equiva- 
lent and their application is dictated by convenience rea- 
sons. The form of kinetic equations ([9|) is symmetric with 
respect to the use of internal and configurational free en- 
ergy. In the right-hand part of Eq. ([9]) the signs are 
chosen basing on solution stability, the internal energy is 
a convex function, and the free energy is a concave one. 
For the well-known Landau-Khalatnikov kinetic equation 



dn dF, 

= -7 



dt 



dn ' 



(10) 



the contradiction is typical. It is written in terms of 
the free energy, but the evolved variable is not an eigen 
argument of the free energy. However, one could show 
that Eq. PU)) can be derived from the second Eq. Q 
in the limit of independence of the defect energy on the 
defect number uyi 0. 

Thus, the Landau-Khalatnikov equation can be only 
generalized to other types of defects with the above con- 
dition satisfied. In the general case, when the defect en- 
ergy depends on the defect number and strongly changes 
in a studied process, for example, in the case of severe 
plastic deformation in metals with grain boundaries as 
defects, it is necessary to use a more general form of the 
kinetic equations Peculiarities of the thermal part 
of the problem can be considered for a special case a 
dislocation generation during the indentation. 
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FIG. 2: Geometry of 2D computer experiment positions of 
atoms at the time 0.21ns. More bright hue corresponds to 
higher potential energy of particles near dislocation cores. 



B. Thermal part of the problem 

Time scans of the wave-like motion of lattice atoms 
can be examined as some signals. Every atom partici- 
pates simultaneously in both the equilibrium and non- 
equilibrium motions. Using the differences in frequency 
properties they can be separated by means of an ordinary 
filtration. 

At present, however, there are no experimental meth- 
ods to registrate the motions of separate atoms. At the 
same time, the molecular dynamics methods allow to si- 
multaneously simulate motions of all atoms of the lattice 
in detail. As an example let us consider the simulation 
of indentation into a 2D copper crystallite consisting of 
80*70 atoms. The indenter (three atoms under arrow) 
moves with a constant velocity of 5m/s (see Fig. [5]). The 
lowest row of atoms is motionless. All lateral faces are 
free. The sample consists of copper atoms interacting 
through the Lennard- Jones potential. A detailed general 
formulation of the problem can be found in Ref. l39l . 

At some moments the dislocations emerge in the 
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FIG. 3: Time scanning: 1-3 total internal, potential and 
kinetic energy, consequently; 4, 5 X- and Y-components of 
velocity of an arbitrary particle in the centre of the model 
(light lines are low-frequency records of particle velocity after 
filtration) 



monatomic specimen, as it is shown in fig. [51 This is 
pointed by bends of the curve 1 in Fig. [3] for time depen- 
dent total internal energy or by the jumps of total poten- 
tial and kinetic energies (curve 2 and 3). The external 
work performed during the indentation increases total 
energies of all kinds. Further, for t = 0.055ns; 0.14ns; 
0.24ns; 0.39ns and so on the potential energy decreases 
sharply. A part of this energy transforms into defect en- 
ergy and the remaining part must dissipate into heat. 

Bottom subplots in Fig. [3] present the time depen- 
dences of the X- and Y-components of velocity for an ar- 
bitrarily chosen particle of the system (near the center of 
the model), which actually are records of its thermal mo- 
tion (incurves 4 and 6). From a formal point of view, the 
thermal motion can be treated as a high-frequency signal 
which is close to a harmonic one for separate interatomic 
(bond) vibrations. The spectrum severely oscillates be- 
cause of the interference of a large set of high-frequency 
phonons (see region 1 in Fig. [4^). In the low- frequency 
region a large peak 2 stands up sharply against the rest 
of the spectrum. One can assume that its presence is 
connected with the regular low-frequency vibrations in 
the system. 

Having done the low-frequency filtration with a run- 
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FIG. 4: Spectrum of complete record for velocity of a particle: 
1 main high-frequency region, 2 low-frequency peak. 



ning average method over 200 time steps one can separate 
the low-frequency vibration modes (bright curves 5 and 
7 against main black signals). To make sure that these 
low-frequency vibrations are connected namely with the 
peak 2 of the spectrum of Fig. 0^, the filtrated record was 
subtracted from the total record and for this difference 
signal the spectrum was calculated again (Fig. HJd). As 
seen, the low-frequency peak vanishes from the spectrum 
of the difference record. 

Comparing curves 2, 3, 5 and 7 one can immediately 
see that the initiation of the low-frequency vibrational 
excitations coincides in time with the generation of dis- 
locations and they can be in a sense associated with 
strongly nonlinear acoustic emission. On the other hand, 
they can be treated as nonequilibrium phonons with am- 
plitudes far exceeding the amplitudes of corresponding 
low-frequency range of the spectrum which would be re- 
alized in thermal equilibrium (compare spectra of fig. 2^ 
and fig. IDs). Thus, using the low- frequency filtration, 
it is possible to separate, in the pure state, the nonequi- 
librium subsystem (light curves 5 and 7) and the equi- 
librium subsystem (the difference signal) in a common 
thermal motion. The acoustic waves inelastically scatter 
by equilibrium high-frequency vibrations and, as a result, 
became damped. Decay of these waves is just a process 
of the nonequilibrium state relaxation. As a result, the 
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FIG. 5: Time scanning: 1 temperature of equilibrium sub- 
system ; 2 temperature of non-equilibrium subsystem. For 
comparison, the scale of curve 2 (the right scale) is 10 times 
magnified. 



energy of the acoustic waves gradually transforms into 
the energy of the equilibrium thermal motion leading to 
a slow increase of equilibrium temperature and entropy. 

During the indentation dislocations emerge repeatedly 
in time. That's why some fraction of the low-frequency 
nonequilibrium wave packets is always present in the in- 
tegral thermal motion. On one hand, they are generated 
as a consequence of the mentioned processes of formation 
of defects, on the other hand, they are constantly going 
to the equilibrium subsystem. By averaging the square 
of particle velocity, that is the kinetic energy, for the 
difference record over time window containing 1000 time 
steps, we obtain an equilibrium temperature-like variable 
(curve 1, Fig. [5]). Such averaging window contains nearly 
20 periods of the high-frequency vibrations, and it seems 
to be enough to smooth of random fluctuations of the 
thermal motion. As seen, the equilibrium temperature 
increases as the energy of nonequilibrium low-frequency 
vibration transforms into the equilibrium high-frequency 
motion. 

As the low-frequency acoustic wave differs from the 
high-frequency phonons only in the time scale then ap- 
plying the same averaging procedure but over a larger 
time window containing in the given example 10000 time 
steps it is possible to definite an analogous variable for 
the non-equilibrium subsystem (curve 2, Fig. [5|). It can 
be given a sense of a nonequilibrium temperature or tem- 
perature of the nonequilibrium subsystem. Thus, the 
digital filtration allows us to separate equilibrium and 
non-equilibrium components in the initial data and to 
calculate their temperature parameters. 

Now, equilibrium S and nonequilibrium S entropies 
can be defined as variables thermodynamically conjugate 
to the two kinds of temperatures described above. Pro- 
duction of nonequilibrium entropy will be described by 
kinetic equation of a form close to ([5]): 



'di 



Here Tq is the temperature of internal sources at dis- 
location generation such that 75X0 is the power of en- 
tropy sources, Tgt is the temperature of internal sinks of 
nonequilibrium heat dissipation such that "t§Tst is the 
power of entropy sinks. It also is a stationary temper- 
ature value, to which the system tends during its relax- 
ation. The last term in Eq. (jlip is a relaxation one, it 
describes transition of the entropy from non-equilibrium 
subsystem to equilibrium one because of the absorption 
of the low-frequency vibrations during their interaction 
with the equilibrium high-frequency vibrations. 



III. THERMODYNAMIC IDENTITY 

The results obtained in the previous section and their 
interpretation allow us to make some useful general- 
izations. In the thermodynamics and theory of self- 
organization, the subdivision of entropy change is well 
known [40] : 



A5 = AeS + AiS, 



(12) 



where A^S is the entropy change owing to heat flux from 
external sources (thermostats) . The heat change is glob- 
ally non-equilibrium but it is a locally equilibrium pro- 
cess. This component is positive at heating and neg- 
ative at cooling. The component AiS is the entropy 
production owing to irreversible phenomena from inter- 
nal sources. Namely, this part of the entropy specifies 
the second law of thermodynamics and is always non- 
negative. The component is of a dual nature. On the 
one hand, this contribution to entropy is caused by non- 
equilibrium and irreversible nature of phenomena. On 
the other hand, this entropy part converts constantly to 
equilibrium form. Therefore, it is pertinent to present 
this part of the entropy as a sum: 



A,S = A^S + Af5, 



(13) 



where AfS is the part of produced entropy succeeded in 
transforming to the equilibrium form during external ac- 
tion, A"^ is the part of produced entropy, which remains 
in the nonequilibrium form and should relax later. The 
former part AfS became indistinguishable from another 
part of the equilibrium entropy (thermal degeneration) 
which had already been existing in the system before the 
action. Thus it makes sense to combine them denoting 



AegS = AeS + AtS, 
so that the total entropy change is: 



AS^ A,gS - 



(14) 



(15) 



-7s (To - Tst -TiS). 



(11) 



In such a form the entropy production has been firstly 
presented in this paper. Usually, one neglects the non- 
equilibrium part of the entropy A^S. Entropy produc- 
tion is accompanied by heat production, which in the first 
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approximation can be presented as a linear combination 
of the entropy change of both kinds: 



A Q = T Ae, ^ + T AI' 5, 



(16) 



where T and T are some coefficients of proportionality 
with a dimension of temperature. Coefficient T has the 
physical sense of the equilibrium temperature. Coeffi- 
cient T can be assigned a non-equilibrium temperature 
sense or a sense of the temperature of non-equilibrium 
subsystem. 

The total change of the internal energy under the exter- 
nal actions and owing to the change of internal structure 
as a consequence of relaxation processes can be written as 
a linear combination of changes of all independent vari- 
ables: 



dlJ — Vaijdelj 



TdS + f l\S + ^ifi /\Hi, (17) 



where more compact notation dS = l^eqS and A5 = 
A"5' is used. The first term is the change of the internal 
energy coinciding with the change of the elastic energy, 
aij is the stress tensor, ef^ is elastic deformation tensor 
(Coleman at all, 1967). The last term stands for the 
change of the internal energy owing to defect subsystems. 
Hi is the number of ^type defects, ipi is the factor having 
the sense of defect energy. Here it is taken into account 
that the internal energy, the elastic deformation and the 
equilibrium entropy are functions of state, therefore their 
decrements are perfect differentials. 

If, following the Landau's idea, one treats the non- 
equilibrium variables T and (pi as generalized thermody- 
namic forces, for which relationships of the type ([3]) 

are true, then in expression (jl7|) . the increments of the 
independent variables AS* and AHi can be replaced by 
their exact differentials: 



JV 



dU = V^o-yde^^ -t- TdS + fdS + ^ ipidHi. 



(19) 



1=1 



Thus, the increment of the internal energy can be ex- 
pressed in a perfect differential form and consequently, we 
can introduce an extended non-equilibrium state speci- 
fied by a set of independent variables e'^j, S, S and Hi. 
A part of these variables efj and S are equilibrated, the 
other part S and Hi are nonequilibrated. In this case, 
the internal energy is a function of these variables 

U^U{elj,S,S,Hi) (20) 

To describe the isothermal phenomena, it is convenient 
to convert to the free energy F = U — TS, for which the 
differential form 



dF = Voij delj 



SdT + TdS - 



N 

E 

1=1 



ipidHi 



(21) 



IV. DISCUSSION 

One can treat relationship (19) as a combination of the 
first and second laws of thermodynamics written in the 
form of identity. Previous attempts to formulate this law 
in the form of an identity are known [4l|, H^l . For exam- 
ple, in the Van's work the Gibbs relationship in the con- 
text of generalized thermodynamics with internal-state 
variables is presented in the form: 



du — Tds 1- ■ de ~ A- da 



(22) 



where u and s are densities of internal energy and en- 
tropy, T is the absolute temperature, e and t are ten- 
sors of deformation (total) and stress, Ais & vector ther- 
modynamically conjugate to the vector of internal pa- 
rameters of. The first two terms are useful attributes of 
the combined first and second laws of thermodynamics; 
the last term A- da > quenches an excess energy. If 
one drops it, one obtains a common expression for the 
law in the form of inequality du < Tds + t- ■ di. In a 
case of the structureless body this inequality arises as a 
common part of the internal energy that has been taken 
in account twice, firstly in the form of irreversible work 
AAn ~ AA — i- ■ de"^ and secondly in the form of inter- 
nal heat released T Af S + f S. And this redundant 
energy is subtracted in the right-hand side of (|22p. 

However, as the total deformation e is irreversible, 
the elastic stresses t can't be calculated as a derivative 
t — du/de. Only the reversible (elastic) part of the de- 
formation can suit the case i = du/de'^ Therefore, 
the structure of relationship (19) more correctly reflects 
the physical nature of the mutual energy transformations 
in internal processes. In relationship (19) two aspects of 
the law of conservation of energy are taken into account 
simultaneously. It is the law of conservation for a given 
reference volume interacting with the environment (the 
1st law of thermodynamics) and the law of conservation 
or the law of transformation of energy in internal degrees 
of freedom AA„ = TAf Af . Moreover, in relation- 
ship (19) both structural aspect (fidHi and dynamic one 
TdS are taken in account too. 

The presentation of the combined first and second laws 
of thermodynamics for a specific case of unstructured 
solids and for isothermal processes in form of Jarzynski's 
identity [4l,|4i|: 



AF W 



kT 



kT' 



(23) 



is also known. Here the angular brackets denote averag- 
ing over all states of the system, AF is the free energy 
increment in an isothermal process, W is the work per- 
formed for the system. If the increment AF and the work 
W are small then, expanding the exponents in a series 
and restricting ourselves to the first terms of the series, 
we obtain: 



is true. 



dF = {W)- 



2kT 



Vaijdsij- 



2kT 



Va,,det,. (24) 
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The first term (W) presents the total work performed 
for the system; the second term {W'^)/2kT presents en- 
ergy dissipated as a result of irreversible internal pro- 
cesses, consequently their difference is the reversible elas- 
tic energy. 

For an isothermal process dT = with neglecting both 
transitional thermal processes dS = and contribution 
from structure of a solid dHi = 0, we obtain from (I^T]) 
dF = VcFijdefj that coincides with the Jarzynski's iden- 
tity in the limit 



V. OUTLINE FOR APPLICATION TO SEVERE 
PLASTIC DEFORMATION 

It is of interest to consider how the obtained above 
general relationships can be applied to solving practi- 
cal problems. One of such problems is a change of the 
internal structure of metals processed by severe plastic 
deformation (SPD) [Ulii, [tj. During SPD the grain 
boundaries are intensively produced resulting in refine- 
ment and fragmentation of the grains. The refinement 
of grains at the expense of multiplication of their bound- 
aries is the main process of SPD and the grain boundaries 
are the main structural defect [i^]. At the same time, the 
main building material for the grain boundaries are dis- 
locations and dislocation piles-up. That's why, for the 
description of SPD these two kinds of defects or the two 
levels of defectiveness need to be taken in account as a 
minimum. With quadratic terms remained in the expan- 
sion of the internal energy 

U = Uo + ipoH- ^ipiH' + ipoH - ^ipiH^ + (bHH, (25) 

we write down the corresponding kinetic equations ([9]) 
for this case: 

OH 



at 

dH 



(26) 



Here variables relating to the dislocations are marked 
by tilde. The last terms describe the mutual infiuence of 
different levels. In an adiabatic limit <^ th the equi- 
librium value of dislocation number (or density) is de- 
rived from stationary condition for the second Eq. pT|) : 



(27) 



This stationary value can be used to exclude the disloca- 
tion number from the first Eq. (pS)) : 

fiJ-f rh 

TH^ = ^n-ViH+^{^o~m^vl^ -ViH. (28) 
at ipi 

Thus, the account of dislocations in the adiabatic ap- 
proach has brought only the renormalization of theory 



constants. Let us consider the next level of approxima- 
tion, we suppose, namely, that the right-hand part of the 
second Eq. ([26|) is not strictly zero, but equals a small 
constant , which depends parametrically on H (quasi- 
adiabatic approximation) : 



1^0 - '^iH + (t)H — e 
Integrating the second Eq. we obtain 



'H 



'Pi 



(29) 



(30) 



Substituting this relationship into the first Eq. (|26l) and 
taking e = Sq + SiH we find 



dH 



at 



= 'Pq -lpYH+ f-eoi - ^eiHt. 



(31) 



' H 



Supposing formally h — v /k and t = e/V, and also ip^ 



ef 



ei — —KTHTfjV'^ /(/) we ge t an equation close to the Ry- 
bin's kinetic equation |48j : 



1 dv 

K de 



—Co + cqi^ + e — i^e. 



(32) 



Note that the first two terms in Eqs. ([3T|) and ([32|) have 
opposite signs. The first term in can be treated as 
annihilation of grain boundaries, going by itself at a con- 
stant speed. The second term describes the generation of 
boundaries, the speed of generation being the higher the 
more boundaries (or fragments) are by this time present 
in a solid. For eo ~ 0.2 [i^, it can be argued that the 
contribution of the first two terms to kinetics of fragmen- 
tation is small. Within this accuracy Eqs. PT|) and ([321) 
coincide. 

Thus, in the framework of the approach, the well- 
known and well-approbated kinetic equation for the de- 
scription of the grain fragmentation during SPD has been 
derived. At the same time, the approach has a more gen- 
eral context that permits to take the influence of thermal 
processes into account and to develop two-mode approx- 
imations with the expansion of thermodynamic poten- 
tials into a higher power series in terms of defectiveness 
of solids (H. 



VI. CONCLUDING REMARKS 

In the paper, a phenomenological approach based on 
generalization of Landau technique to cases of quick- 
passing nonequilibrium phenomena, for example, severe 
plastic deformation, is considered. For quick-passing pro- 
cesses thermal fluctuations have no time to exert essen- 
tial influence and it is possible to consider the problem 
in the mean-field approximation. In the centre of the 
approach there is not an abstract order parameter but a 
structural defect with its real parameters the quantity 
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(density) of defects and the average energy of a defect. 
The analysis from these positions on the example of a 
version of the theory of vacancies shows that the known 
Landau-Khalatnikov kinetic equation can be used in this 
case only if the dependence of vacancy energy on den- 
sity of vacancies is neglected. In the case of an arbitrary 
dependence of vacancy energy on density of vacancies it 
is necessary to use a more general form of kinetic equa- 
tions symmetric with respect to using the internal and 
free energies. In this case, the density of defects and 
defect energy are related by a symmetric differential de- 
pendence of type ([5]) and As the defect energy in the 
steady state is not equal to zero then the extreme prin- 
ciple of equality to zero of the derivative of free energy 
with respect to order parameter (in our case the num- 
ber or the density of defects) breaks down. This prin- 
ciple need be substituted with principle of the tendency 
to a steady state. Steady-state characteristics can not 
be determined in the framework of phenomenological ap- 
proach, here statistical and microscopic approaches are 
required. 

A form of kinetic equations can be generalized to all 
types of regularly or randomly distributed defects and to 
all types of nonequilibrium parameters including thermal 
variables. On an example of computer experiment it has 
been shown that in the case of quick-passing processes the 
irreversibility is not so much connected with fluctuation 
phenomena, as with dynamical ones of acoustic emission 
type at the formation and motion of defects. In this 
case, the acoustic emission over all structural levels is a 
part of the relaxation process [H, IH, [s^ that permits to 
describe the entropy generation (|lip in the same form as 



generation of defects ([9]). 

The obtained results are generalized in the form of a 
thermodynamic identity which properly combines the 1st 
and the 2nd laws of thermodynamics in a finite-difference 
form. The known Jarzynski thermodynamic identity for 
isothermal processes expressed in terms of free energy 
coincides, in going to finite differences, with the proposed 
identity. 

To illustrate the applicability of the approach a prob- 
lem of description of nonequilibrium phenomena occur- 
ring under severe plastic deformation is considered. In 
the framework of phenomenological one-mode (approxi- 
mation quadratic with respect to energy) model taking 
into account the two levels of defect structure dislo- 
cations and grain boundaries, the known for this case 
kinetic Rybin equation is derived [i^. The approach 
permits to consider more complex dependences of ther- 
modynamic potentials on basic parameters of the model. 
For example, in the framework of two-mode approxima- 
tion using the polynomials of the fourth power for pre- 
sentation of thermodynamic potentials the severe plas- 
tic deformation can be represented as a structural phase 
transition of the 1st type 
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